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Abstract . This is an expanded version of the talk given at the conference "Constructive Functions Tech-04" . We survey 
some recent results on canonical representation and asymptotic behavior of polynomials orthogonal on the unit circle with 
respect to an analytic weight. These results are obtained using the steepest descent method based on the Riemann-Hilbert 
characterization of these polynomials. 
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1. Introduction. During the Fall Semester of 2003 I was visiting the Department of Mathematics of 
the Vanderbilt University, where I had the opportunity to continue my collaboration with Ed Saff. I was 
very excited with the evolution of the Riemann-Hilbert approach to the asymptotic analysis of orthogonal 
polynomials, and discussed extensively with Ed the new perspectives. He was the one who posed the question: 
can this method tell us anything new about such a classical object as the orthogonal polynomials on the unit 
circle (OPUC, known also as Szego polynomials), in particular, about their zeros? The question was more 
on a skeptical side. I was aware of some previous work of the founders of the method, p^, [3], but none of 
these papers was focused on the description of the zeros of the OPUC's. So, we started to work and realized 
that we were able to find curious facts even in the simplest situations. Later on I visited Ken McLaughlin, 
at that time in Chapel Hill. A two-day discussion at Strong Cafe (a recommended place) was crucial, and 
Ken joined the team. This paper is a short and informal report on some of the advances we have had so far. 

Let me introduce some notation and describe the setting. For r > 0, denote H) r = {z € C : \z\ < r}, 
and T r = {z £ C : \z\ = r}. A positive measure /i on Ti has the Lebesgue-Radon-Nikodym decomposition 

d/i(z) = w(z) \dz\ + d/i s (1-1) 

where fj, s is the singular part of fj, with respect to the Lebesgue measure on Ti. Throughout, we will consider 
measures satisfying the Szego condition 

logw(z) |g?z| > — oo , 

Ti 

allowing to define the Szego function (see e.g. ^3 Ch. X, §10.2]): 

1 f 2 \ ,^e id +z 



D(w;z)=exp[— \ og w{e w ) -^—dO ) . (1.2) 



This function is piecewise analytic and non- vanishing, defined for \z\ ^ 1, and we will denote by D\ and D c 
its values for \z\ < 1 and \z\ > 1, respectively, given by formula 11.2J1 . It is easy to verify that 



z J D c (w; z) ' 

and for the boundary values we have 



A ( w;l ) = _ * . , \z\>l, (1.3) 



, , D\(w:z) 1 ™ / \ 

= n \ 4 = T7T7 VRF' z e T i ■ L4 

D c {w]z) \D e (w;z)\ 2 
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The first equality in l|1.4(l can be regarded as a Wiener-Hopf factorization of the weight w, which is a key 
fact for the forthcoming analysis. 

For a nontrivial positive measure fj, on Ti there exists a unique sequence of polynomials (p n (z) = n n z n + 
lower degree terms, K n > 0, such that 



ip n (z)(f m {z) dfi(z) = S mn , m,n = 0,1,... (1.5) 

We denote by & n (z) = tpn{z)/n n the corresponding monic orthogonal polynomials. They satisfy the Szego 
recurrence 

$„+i(z) = z$ n (z) -a^$* n {z) , $0(2) = 1, 

where we use the standard notation $* (z) = z n <& n (l/z). The parameters a n = — $„+i(0) are called Verblun- 
sky coefficients (also reflection coefficients or Schur parameters) and satisfy a n G Di for n = 0, 1, 2, 

Furthermore, /1 <-> {a„} is a bijection; the map /1 — > {a„} is an inverse problem, and it is known to be 
difficult. In JHj there is a thorough discussion of several techniques to tackle this problem. Recently, the 
Riemann-Hilbert approach, not described in JH|, proved to be very promising in this context also (see pQ, 
0, El)- The main goal of this paper is a further discussion of how this method can shed new light on 
the study of the asymptotics of the Szego polynomials. We are not going to provide detailed proofs that 
can be found elsewhere (the references are included), the aim is to show the method in action and to discuss 
some new results in two apparently simple situations. 

The structure of the paper is as follows. Sectional is devoted to the case when D(w; •) is non vanishing 
and has an analytic extension across Ti. The main role here is played by the scattering function 1 

S(w; z) = Di(w; z)D c (w; z) , (1-6) 

meromorphic in an annulus, containing Ti, and which, via its iterated Cauchy transforms, allows to write 
some canonical series representing y n 's. In this situation convergence is always exponentially fast, and the 
Riemann-Hilbert analysis is particularly simple and transparent. Only some of the multiple corollaries of the 
canonical representation are discussed; a more thorough analysis is contained in |12| . In Section [3] we look 
at the situation when the original analytic and nonvanishing weight has been modified by a factor having 
a finite number of zeros on the unit circle. Now the behavior of the zeros of 'I'n's is qualitatively different: 
most of them cluster at Ti, and only a finite number stays within the disc 13%. The method of Section 
must be modified now in order to handle the zeros of the weight: a local analysis plays the major role. 
The exposition here is much more sketchy; in this sense, more than a detailed view this window gives us a 
glimpse of the possible techniques and results. At any rate, the main goal is to persuade the reader that the 
Riemann-Hilbert analysis is a powerfull technique, that deserves to be in the toolbox of those interested in 
orthogonal polynomials on the unit circle. 

2. Analytic and nonvanishing weight. Any analysis of OPUC's can be started either from the 
orthogonality measure /i or from the sequence of the Verblunsky coefficients {a„}. Let us assume that the 
sequence of the Verblunsky coefficients has an exponential decay: 

p= Em |a„| 1/n . (2.1) 

n — »oo 

Nevai-Totik |15| proved that this situation is characterized by the following conditions on fi: in the decom- 
position fi s = 0, measure fj, satisfies the Szego condition, and 

p = inf{0 < r < 1 : D c (w; z) is holomorphic in \z\ > r} . (2-2) 

1 Function 1/S is denoted in 1171 by r, and in 1181 Section 6.2] by b. It corresponds also to the scattering matrix in 0. I 
prefer to follow the notation of 1121 . 
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Taking into account ljl.3|l we see that the first identity in (|1.4|l can be regarded as an analytic extension of 
the weight w. With this definition of w (which we use in the sequel) we can say equivalently that 

p = inf{0 < r < 1 : l/w(z) is holomorphic in r < \z\ < 1/r} , (2-3) 

and both circles T p and Ti/ p contain singularities of l/w. In this situation the well-known Szego asymptotic 
formula 

Urn <p* n (z) = —-r- — r, ZGDi, 
n— >oo L)\\W\Z) 

can be continued analytically through the unit circle Ti and is valid locally uniformly in H\/ p . It shows that 
the number of zeros of {^p n } on compact subsets of EDi \ B p remains uniformly bounded, and these zeros 
are attracted by the zeros of D e in p < \z\ < 1 (Nevai-Totik points in the terminology of B. Simon |16|h 
Numerical experiments show that the vast majority of zeros gather at the "critical circle" T p . This fact was 
justified theoretically by Mhaskar and Saff ^1], who using potential theory arguments showed that for any 
subsequence {n^} C N satisfying 

p = lim|a„J 1 /"S 

the zeros of {(p„ k+ i} distribute asymptotically uniformly in the weak-star sense on T p . 

The behavior of the zeros inside D p can be intriguing. Although approaching in mass the critical circle 
T p as predicted by Mhaskar and Saff, some of them still may remain inside and follow interesting patterns. 
Even the convergence to the circle T p is different for different measures. Can we give a full description of this 
behavior in terms of the weight of orthogonality? The answer is positive, and the description will involve 
a sequence of iterates of certain Hankel and Toeplitz operators with symbols depending on the scattering 
function S introduced in ljl.6|l . It is a consequence of a canonical representation of the Szego polynomials, 
found by means of the Riemann-Hilbert characterization. 

2.1. Steepest descent analysis and canonical representation for orthogonal polynomials. 

The starting point of all the analysis is the fact that under assumption (|2.1(l conditions (|1.5|l can be rewritten 
in terms of a non-hermitian orthogonality for ip n and 95* : 

(p n {z)z n - k - 1 ^-dz = 0, for k = 0,1,.. 1, 

z n 

k w(z)^_jo, fc = 0,l,...,n-2, 



Ti z 



if Kn-i, k — n — 1 . 



Here and in what follows, all the circles T a , a > 0, are oriented counterclockwise; with this orientation we 
talk about the "+" and the "— " side of T Q referring to its inner and outer boundary points, respectively. 
Analogously, / + and /_ are the corresponding boundary values on T Q for any function / for which these limits 
exist. By standard arguments (see e.g. PP or [7], as well as the seminal paper where the Riemann-Hilbert 
approach to orthogonal polynomials started), 

/ j_ r * n (tMt)dt \ 

W 2™ J Tl t»{t-z) 

is a unique solution of the following Riemann-Hilbert problem: Y is holomorphic in C \ Ti, 

y+(t) = F_(t) (J w(< j A ") , Z£T U and \imY(z)( Z Q l °„) = / , (2.4) 
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where I is the 2x2 identity matrix. 

This is the starting position for the steepest descent analysis as described in [J] (see also JD]), which 
consists in performing a series of explicit and reversible steps in order to arrive at an equivalent problem, 
which is solvable, at least in an asymptotic sense. Since these steps are almost standard, they will be 

described very schematically. We will use the following notation: 03 ' 







-1 



is the Pauli matrix, and 



for any non-zero x and integer m, x a3 
Step 1: Define 



x 



and x m ° 3 = {x rn ) c 



H(z) 





„ \z-™\ ifM>l, 



if M < 1, 



(2.5) 



and put T(z) = Y(z) H(z). Then T is holomorphic in C \ Ti; this transformation normalizes the behavior 
at infinity: lirn^oo T(z) = I. The price we pay is the oscillatory behavior of the new jump matrix on Ti: 



T + (t)=T_(i) 



t n w(t) 
t~ n 



We get rid of these oscillations in the next transformation, taking advantage of the analyticity of its entries 
in the annulus. 

Step 2: Choose an arbitrary r, p < r < 1, that we fix for what follows; it determines the regions (Figure 

IP) 

fi = {z : \z\ < r}, fl OD = {z: \z\ > 1/r}, 

0+ = {z : r < \z\ < 1}, = {z : 1< \z\ < 1/r}. 

Define U(z) = T(z)K(z), where 




Fig. 2.1. Opening lenses. 
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K{z) = i 



(I, 

1 °V' 

z n /w(z) 1 J 
1 N 



, if z e f2 + , 

, ifzefi_. 



(2.6) 



Then {/ is holomorphic in C \ (T r U Ti U T\/r)- We have not modified its behavior at infinity, but now 

U+{t) = C7_(f) Ju(t), t e (T P U Ti U T 1/r ) , 



where 



Ju(t) = < 



w(t) 
-l/w(t) / 

1 0\ 
t n /w{t) 1 j ' 

1 0\ 
l/(t n w(t)) 1 



, ifiGTi, 
if t € T r , 
if i € T 1/r 



The jump on T r and Ti/ r is exponentially close to the identity, which is convenient to our purposes. We 
have to deal now with the relevant jump on the unit circle. 

Step 3: The Szego functions D\ and D c have been introduced in 11.2J1 . Define 



1 / 1 f 2 * 

— — — = D e (w; oo) = exp -— / log 

A(w;0) V 4?r Jo 



w {e % °)d6 > 0. 



Hence, if we introduce the geometric mean 



Q[w] = exp (J- J 2 J log (to (e w ) dB) 



(2.7) 



(2.8) 



then t — {Q\w\) 1 ^ 2 . 

It is straightforward to check that the piece-wise analytic matrix-valued function 



N(z) = N(w; z) = < 



D c {w-z)/t \ 

t/D c {w-z)) 
A(w;z)/r N 
-t/Di(w;z) 



if \z\ > 1, 
if \z\ < 1, 



(2.9) 



is invertible, and has the same jumps on Ti as U(z). This motivates to make a new transformation, defining 
S(z) = U(z)N- 1 (z). Matrix S is holomorphic in C \ (T r U Ti U T 1/r ), 



lim S(z) = I, 



and 



S+(t) = S-(t)J s (t), t e T r U Ti U T 1/r 



(2.10) 
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where 



J s (t) = N-JuN- 1 = < 



1 -t n S{w;t)/T 2 
1 

1 0' 
\T 2 /(t n S(w;t)) 1 



if t S Ti, 
if t e T r , 

if t e Ti 



(2-11) 



A 



Our main character, 5, has entered the picture! 
Summarizing, we have 



Y(z) = S{z)N(z)K- 1 (z)H- 1 {z) . 



(2.12) 



Here H, K and N are explicitly defined in (|2.5p . i|2.6f) and (|2.9I) . respectively. About 5 we know only that 
it is piece-wise analytic and satisfies the jump condition H2.10fl - 12.llfl . A feature of this situation is that we 
can write a formula for S in terms of a series of iterates of some Cauchy operators acting on the space of 
holomorphic functions in C \ (T r U T]/ r ) with continuous boundary values. Indeed, let us denote 



S 



Su S12 
5*21 S 22 



and look at the equations given by the jumps on T r . According to (|2.11fi . the first column is analytic across 
this circle, while the second column has an additive jump equal to the first column times — t n S(w; t)/T 2 . So, 
if we define the operator 



M l n (f)(z) 



2ni t 2 



where /_ denotes the exterior boundary values of the function / on T r , then taking into account the behavior 
at infinity and Sokhotsky-Plcmelj's theorem we get that 



Si 



M^Sn), S22 = l + M l n (S 2 i) 



(2.13) 



Analogously, with 



i 



S{w;t) t n (t-z) 



dt, 



where /_ denotes the exterior boundary values of the function / on Ti/ r , we have again 

Su = l + M e n (Si 2 ), S 2 i=M c n {S 22 ). (2.14) 

By If2.13fl . (j2.14fl . functions Sij satisfy the following integral equations: 

(I - Ml o M^Sii = 1 , (I - M„ a Ml)S 12 = M [ n {\) , 
(I - M e n c M [ n )S 21 = M c n (l) , (/ - M\ o M c n )S 22 = 1 , 

where / is the identity operator. Straightforward bounds show that there exists a constant C > depending 
on r only, such that 



\Ml(f)(z)\ < Cr 
\Ml(f){z)\<Cr 



n Wf-h 



Tr 



z\-r\- 

n \\f-hi /r 

\\z\-l/r\ 



(2.15) 



T 1/r 
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where || • \\j is the sup-norm on 7; in the sequel we use C to denote some irrelevant constants, different in 
each appearance, whose dependence or independence on the parameters will be stated explicitly. Thus, we 
can invert these operators using convergent Neumann series, 



Ski = ( ]T {Ml o M\) k (1) , Ska = (jr (Ml o M c n f o A< (1) , 

\fe=0 / \fe=0 / 

S 21 =|X;KoK) fc oM c J (1), S 22 = (f^(M) l oMl) k ) (1). 



\fc=0 / \fc=0 



In order to make this somewhat more explicit, let us define recursively two sequence of functions: 

/r = l, fP^MUl), and fW^Ml(fi 2k -V) , /^ +1 > d = f MUil 2fe) ) , fc^N, 

ff (°) d = f l, 9 W d = f M^(l), and ff^^KCffg 1 *- 1 )), 9 r +1)d = f M e n (^), keE. 

Then 

OO OO 

Sii(n^) = ^4 2fe )( Z ), S 21 (n;z) = ]T 5 i 2fe+1 >(z), for |z| ^ 1/r, 

fc=0 fc=0 
00 00 

S 12 (n;z) = Y,fn k+1 Hz), S 22 {n-z) = Y.9 { n k \z), for |«| + r . 

k=0 k=0 

All these series are uniformly and absolutely convergent in their domains. Moreover, straightforward bounds 
on the Cauchy transform show that there exists a constant C > depending on r only, such that for i = 1,2, 

C C 
\Sn(z)\ < _ , for \z\ y£ l/r, and \S&(z)\ < ,, , _ , , for \z\^r. 

We emphasize that in each of the regions above functions Sij have their own meaning, and are not obtained 
in general by analytic continuation from one domain to another. 

Once we have computed S, we may replace its expression in (|2.12() in order to find Y. It should be 
performed independently in each region. For instance, in the domain Cloo, N(z) = (D e (w; z)/r) <T3 , K(z) = /, 
H(z) = z- na3 , so that 



Y(z) = S(z)(z n D e {w;z)/r) 



ff3 



Sn(n; z)z n D c (w; z)/t Si 2 (n; z)T/(z n D c (w; z)) 
S 2 i(n; z)z n D c (w; z)/t S 2 2(n; z)r/(z n D c (w; z)) 



Analogous computations are easily completed in the rest of the regions. 

All the information about the parameters of the orthogonal polynomials is codified in the first column 
of Y: its (1, 1) entry gives us <!>„, that evaluated at z = yields the Verblunsky coefficients, while the (2, 1) 
entry at z = is related to the leading coefficient k: Y2i(0) = — 2ttk^ 1 _ 1 . Hence, we have obtained the 
following theorem, that has been proved in |12j : 

Theorem 2.1. Let w be a strictly positive analytic weight on the unit circle T\, the constant p as defined 
in (|2.2|) - ||2.3|l , and constant r with p < r < 1 fixed. Then with the notations introduced above, for every 
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n £ N sufficiently large the following formulas hold: 

r -1 z n D c (w; z) Sn{n\ z) 

l) <f> n (z) = < 



if\z\ > 1/r; 

T- 1 z n D e (w; z) Su(n; z) - T ^f n ' , z / r < | z | < i/ r; 



r5i 2 (n;z) 

m) 5^" = t 2 Si 2 (7J+ 1;0). 
m) k 2 = — S 2 2(n + 1;0). 

Z7T 



A(w;z) 



«/ |z| < r . 



(2.16) 

(2.17) 
(2.18) 



Remark: It is easily seen that the method we have just described is valid also if we replace the condition 
of positivity of w on Ti by the requirement that its winding number on Ti is zero. In such a case we can 
assure that deg <!>„ = n only for sufficiently large n's, but the rest of the argument remains the same. 

Before we analyze some implications of these formulas let us look more carefully at the scattering function 
S and at the iterates of its Cauchy transforms used in the definition of S. Following notation of |18l Section 
6.2], let 

logw(z) = ^L k z k 

fcez 

be the Laurent expansion (equivalently, the Fourier scries) for logic. Then straightforward computation 
shows that 



Let 



S(w; z) = exp ]T [L k 

\k=l 



+00 



z k -L_ k z- k 



S(w;z)= {S) k z k and 



exp 



1 



Y J [LkZ k -L k z k 



\k=l 



k= — oo 



S(w; z) 



+00 , 

£ (s 



k=—oo 



(2.19) 



be the Laurent expansions of S and 1/5 in the annulus p < \z\ < 1/p, respectively. Taking into account that 
\S(w; z)\ = 1 on Ti, 



+00 



5] =( 5 )-A> and J2 ( S )k +m (S) k = 



k= — oo 



1, if m = 0, 

0, ifm€Z\ {0} . 



(2.20) 



Denote by ffi the Hardy class, and H 2 _ = L 2 H 2 :. Let us denote by V+ and V- the Riesz projections 
onto and H 2 , respectively, and a„ (z) = z n S(w; z). Then 



+00 



V + (a n )(z)= 



k-\-n 



k——n 



k<C—n 



In particular, 



-Hoc- 



1 1 7^ 

--V + (* n )(z) = --5 ]T (S) fe ^+"etf 2 , if|z|< 

k— — n 

^V-{a n )[z) = i ( S )k* k+n € Hi , if 



z > r . 



k<—n 
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The series in the right hand side converge locally uniformly. Observe also that 



/«(0) = -l(<S)_ n = -i(i^. (2.21) 

If we introduce the following composition of Hankcl and Tocplitz operators, having cr^ 1 as a symbol, 
Hi: H 2 _^H 2 ± , given by H+(f) = -V+(a n P- K 1 /)) , = P-M-fcV)) , 

then 



,(2 fc +i) = ^ N<r, 

Wi^), if|*|>r, 



where /„ * represents the values of /^ 2fe 1 in C \ D r . 
For we can obtain analogous formulas: 

we i = /u 1 r = J- T_2:P +K^ 13 )(2) = -p+(o»p+(0)(z) » ^ N < r, 

■ 9 " U n(9n ' \r- 2 V-(a n g^)(z)=V-(a n V + (a-'))(z), if\z\>r, 

and 

12k) 



(2fc+2) = I W+(5„,J), if|z|<r, 
For instance, taking into account (|2.2Q|1 . for \z\ < r, 



^ 2) w = -n(^ + (0)W = -^+( ( E (^^"l (e(|) 

y \j — — oc J \k—n ^ 
+oo / 

= - e W^'-^-E E ^ 

fc<j, fc<— n m— Y^— — 71 

and by identities in (|2.20(1 we can rewrite last formula as 

-+-00 / \ 



!+^)=E E(%J^ z™, \z\<r. (2.22) 



m=0 \fc>— n 



Remarks: There are some further relations and equivalent expressions that an interested reader can 
easily derive. For instance, if we introduce the operator T n on 1? with kernel 



s=0 ^ ' s+n+l+i k<-n-j ^ ' 



i-j-k 



then 



P+T„(/), if|z|<l/r, 
-P-T n (f), if|z|>l/r, 



and we can obtain expressions for fn k ^ and g\i k+1 \ 
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2.2. Asymptotic behavior of OPUC. Representation (|2.16(l is asymptotic in nature. Using lj2.15|l . 
it is immediate to show that for all sufficiently large n and for N = 0, 1, 2, . . . , 



N 



S 11 (n;z)-Y / fL 2k) ^) 



k=0 



< 



c 



\z\-l/r\ 



M^l/r, 



where the constant C depends only on r and N, but neither on n nor on z. Analogously, 

N 



S 12 (n;z)-£/r +1) M 



k=0 



< 



C 



„(2JV+3)n 



(2.23) 



where C has a similar meaning as above. These bounds show that (|2.16J) allows us to obtain approximations 
of {$„} of an arbitrarily high order. We will concentrate only on the most interesting domain including the 
critical circle T p and its interior (for a full analysis, check |12|). 

Let us discuss the consequences of truncating Sn and S\2 in l|2.16|) at their first terms, 



S n (n;z) = 1 + G(r 2 "), S 12 (n; z) = - 



27TIT 2 



S(w;t) V 

t - z 



dt + 0(r in ), 



imposing some additional conditions on the analytic continuation of our weight w (or function S). We 
assume first that the critical circle T p contains only isolated singularities in a finite number. 

Theorem 2.2 ( \1'2\). Assume that there exists < p' < p such that D e can be continued to the exterior 
of the circle T p /, as an analytic function whose only singularities are on the circle T p , and these are all 
isolated. Denote by ax, . . . , a u the singularities (whose number is finite) of D c on T p , 

|oi| = • • • = \a u \ = p. 

Then for p < r' < r there exist constants < 8 = S(r') < 1 and C — C(r') < +oo such that for p' < \z\ < r' 
and n € N, 



®n (z) - z n n c ; \ - ' / es s W ' * 7 

D c {w\ oo) A(w;z) ^t=a fc V i 



< C (p n <5" + r 3n ) . 



(2.24) 



Furthermore, for every compact set K C D p there exist constants < <5 = <5(i ; i') < 1 and C = C(K) < +oo 
smc/i ttai for z € iC, 



A(w;0) 

A(w;^) ^*= Q " fe 



^ res ; 



t-2 



< C (p n S n + r 3n ) 



(2.25) 



If D e can be continued as an analytic function with a finite number of isolated singularities to whole 
disc Di, then we may take 5 = in the right hand sides in (|2.24|l - l|2.25(l . Otherwise the right hand side in 
(|2.25|l may be replaced by an estimate of the form Cp n S n . 

In order to isolate the zeros of $ n , one must be able to analyze the approximation to <&„ afforded by 
(|2.24|) and (|2.25l) . For example, zero- free regions may be determined by (i) establishing zero-free regions for 
the approximation, and (ii) bounding $„ away from zero using the error estimates. Similarly, isolating the 
zeros can be done by first isolating the zeros of the approximation, and then using a Rouche' type argument 
for $„. 

This theorem tells us that in general all the relevant information for the asymptotics of <!>„ 's in D r comes 
from the singularities of the exterior Szego function D e on T p (that is, from the first singularities of D e we 
meet continuing it analytically inside the unit disc), and reduces the asymptotic analysis of $n's (at least, 
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in the first approximation) to the study of the behavior of the corresponding residues. In the case when all 
the singular points that we met on T p are poles, this analysis is more or less straightforward. 

Definition 2.3. Let a £ Di be a pole of a function f analytic in \z\ < 1. We denote by mult z=a /(2;) 
its multiplicity and say that a is a dominant pole of f if for any other singularity b of f , either \a\ > \b\ or 
\a\ = \b\, but then b is also a pole and 

mult f(z) > mult f(z) ■ 

z=a z=b 

In the sequel we use the following notation: for a £ C and e > 0, 

B s (a) = {z£C: \z - a\ < e} . (2.26) 

Theorem 2.4 f |12|). Assume that there exists < p' < p such that D e can be continued to the 
exterior of the circle T p > , as a meromorphic function whose only singularities are on the circle T p . Denote 
by oi, . . . , a u the poles (whose number is finite) of D e on T p , and assume that the dominant poles of D c are 
di, . . . , at, i < U, and their multiplicity is m. 

Let e > 0. Then for p' < \z\ < r - e, z U^ =1 S £ (afe) ; and n e N, 

D c (w;z) ^ n ( Dj{w;0) ( 71 ^ „n-m+i Di{w; a k ) D c (w; a k ) 

k=l 



$ n{z) = J±pj± z n + ^^y2 » a ~-m + l i^Wj +Mz)> (2 . 27) 

L>c(w;oo) Di(w;z) \m- I J a k - z 



where D e (w] ak) = lim z ^ afc D c (w; z)(z — ak) m , k = 1, . . . , I. There exist a constant < C < +oo independent 
of n and e, and a constant < 5 = S(e) < 1, such that 



C (p n S n + r 3n ) , if m = 1 . 

p n , if m > 2 . 



IM*)I<< c ^ m _„„ 



Furthermore, for every compact set K C D p there exists a constant C = C{K) < oo such that for z £ K , 
and n £ N, 



T Di{w,z) f n \ 1 _ A(w;afc)-P c (w;a fc ) 2m (n~m+i)e k 

a n- m +l y m _ 1 J afc _ z 




if m = 1 , 

(2.28) 

— , i/ to > 2 , 
n 



f?i = l, and = (argafc — argai) , k = 2, ...,£. (2.29) 

27T 

In particular, on every compact set K C D p; /or aZ/ sufficiently large n polynomials $ n can /iawe at most 
£ — 1 zeros, counting their multiplicities. 

Observe that this result is applicable to weights of the form w(z) = \R(z)S(z)\ 2 , z £ Ti, where R is a 
rational function with at least one zero on T p (or one pole on T 1 / p ), and S is any function holomorphic and 
^ in any annulus, containing {p < \z\ < 1/p}. 

By means of (I2.27f) we may show also that under assumptions of Theorem 12 .41 the vast majority of zeros 
approaching the critical circle T p does it in an organized way, exhibiting an equidistribution pattern: if zeros 

(n) 

Zj of can be numbered in such a way that, roughly speaking, 

\4 n) \ = p(i + -io S ( \) + o( 1 - 

\ n \m — 1 / \n 
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and 



arg | 



O 



(again, we refer the reader to ^2 for details). This is also an analogue of the interlacing property of the 
zeros of orthogonal polynomials on the real line. Moreover, H2.28(l allows us to describe the accumulation 
set of zeros of {$„}'s inside D p . For instance, if all Ok's in H2.29|l are rational, this set is discrete and 
finite. Otherwise, as it follows from Kronecker-Weyl theorem, it can be a diameter of D p or even fill a two 
dimensional domain. 

The situation gets much more complicated if the first singularity that we meet continuing D c analytically 
inside is more severe. Consider the simplest example of an essential singularity on T p : 



;(*) = 



exp 



1 



p-t 



t e Ti 



(2.30) 



with < p < 1. Observe that its inverse, 1/w, satisfies also the conditions of Theorem 12.21 However, the 
behavior of the zeros of the OPUC for w and 1/w is qualitatively different, check Figure l2~2l 





Fig. 2.2. Zeros of $30 for weights w (left) and 1/w (right), with w given in l2~30l and p = 1/2. 

In a few words, the explanation for this phenomenon is the following. Observe that in the case of an 
essential singularity of S the asymptotic behavior of the Cauchy transform 

Jn v ; 2niT 2 J Tr t-z 

is not as simple as when the only singular points are poles. In fact, the leading term of the asymptotics will 
come now from a dominant saddle point of 

* n (t) = \ogt+- \ogS(w;z) 
n 
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lying close to the singular point of S, t = t+, which is the solution of the equation 

\ - — j— ( — + -i- ) = satisfying t+ = p + + O ( - ) , n oo , (2.31) 

i n + 1 \ pi — 1 t — pj y n + 1 \ n / 

where we take the positive square root. It is possible to show that the zeros of the orthogonal polynomials 
(at least those not staying to close to z = p), will approach the level curve 

1 / 1 p 3 / 4 \ 

Re(vM^) ~ *„(*+)) - - log , (2.32) 

and the error decreases with n. However, for the weight l/w we will have two dominant saddle points, and 
the different structure of the level curve (|2.32|l for weights w and l/w explains the different result of the 
numerical experiments (see Figure l|2.3f) ). 




-1 -0.5 0.5 -1 -0.5 0.5 



Fig. 2.3. Left: zeros o/$„ for w given in 12.301 with p = 1/2 and n = 30, along with the level curves Re(^ n (^)— = 
0, Re(\P n («) — = i log (^T^TTl) and ^(^bW ~ *"(*+)) = 0- Right: the same, but for the weight l/w. Observe 

that the level curve 12.321 has now two components. 



2.3. Verblunsky and leading coefficients. Let us finish this Section with some comments about the 
behavior of the notorious coefficients related to the OPUC. 

Evaluating the polynomial <& n +i or any of its approximations at the origin we obtain information about 
the Verblunsky coefficients. For instance, a combination of H2.17J1 . i|2.21[l . and H2.23(l yields the following 
estimate of the Verblunsky coefficients a n : 

Proposition 2.5 f|12|). Let w be a strictly positive analytic weight on the unit circle T%. With the 
notation introduced above and for each n G N, 

+ 0(r 3n ), (2.33) 

n+l 

where (l/S) n +i is the corresponding Laurent coefficient ofl/S in l|2.19|l . 
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This fact has the following reading: consider the generating function of the Verblunsky coefficients, 

oo 

G(z) = a nZ T > 
n=0 



Then the Maclaurin series of G and V+(— z/S(z)) match up to the 0(r 3n ) term. In consequence, we have 
Proposition 2.6. Function 

G(z) 



S(w; z) ' 

defined in a neighborhood o/Ti, can be continued as a holomorphic function to the annulus 1 < \z\ < 1/p 3 - 
This fact has been established independently by Simon ^7] and Deift and Ostensson 6 . 

If the only singularities on T p are dominant poles, we can use formula (|2.28() in order to derive the 
asymptotic behavior of a„'s: 

Proposition 2.7 ( 12 j). Under assumptions of Theorem \2.4\ the Verblunsky coefficients satisfy 

a n = - 

fc=i 



n+l\ ~ n - m+ln , , , \0{p n S n ), ifm = l 



E , «r +1 AKa fe )5 c K« t )+ ln , 

^ — ' V m — 1/ I [n m 1 p n ) , z/ to > 2 . 



That is, in the situation when the first singularities of D e met during its analytic continuation inside are 
only poles, the Verblunsky coefficients are asymptotically equal to a combination of competing exponential 
functions with coefficients that are polynomials in n. We can compare it with the case of the essential 
singularity considered before: for the weight w given in H2.30|) . 

where t + — > p is given by the equation (j2.31|l . 

A similar analysis can be carried out for the asymptotic expansion of the leading coefficients K n . By 
(l2~m 

^ = ^S 22 («+1;0) = ^ (l+ 9^,(0) +0(r in )) . 

Taking into account l)2.22[l . we arrive at 

Proposition 2.8. For the leading coefficient n n the following formula holds: 

< = ^ E l(S) fc | 2 + 0(r 4 ") = ^ [l- £ \(S) k \A+0(r in ), (2.34) 

fc>— n— 1 V fc<— n— 1 I 

where {S) k 's are the coefficients of the Laurent expansion of S in l|2.19|l . Observe that we can write this 
identity also in terms of the Riesz projections: 

4=r 2 \\P- {°nli)t* {W) +0{r^). 

Formula l|2.34jl shows that 

lim4 = l!, an d a = g \ {S )_ n _^ + 0(r 4 «) , 
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in accordance with Q2.33|l and the well known fact that 



1 



Summarizing, we see that the Laurent coefficients of S (or of 1/5) contain surprisingly good approxi- 
mations of two main parameters of the OPUC: they match asymptotically the Verblunsky coefficients, and 
the partial sums of the squares of their absolute values represent (up to a normalizing constant) the leading 
coefficient of the orthonormal polynomials. 

3. Weight with zeros on TV Let us analyze the change of the behavior of the orthogonal polynomials 
if we allow zeros of the weight on the unit circle. In other words, we consider now a weight of the form 

m 

W(z) d = f w(z) Y[\z- a k \ 2 ^ , z G Ti , (3.1) 
fe=l 

where a k € Ti, (3 k > 0, k — l,...,m, and w is an analytic and positive weight on Ti, such as considered 
in Section Without loss of generality we assume that w is analytic and non- vanishing in the annulus 
p<\z\< l/p. 

According to Nevai and Totik JS], the Verblunsky coefficients no longer have an exponential decay, 
neither the bulk of zeros accumulate on an inner circle, but how many of them stay inside? And for those 
approaching the unit circle, does the rate depend on the "orders" /3fc? And how can we extend the Riemann- 
Hilbert method, that so nicely worked for us in the analytic situation, to the case of a weight of the form 

(EH)? 

3.1. Steepest descent analysis. Matrix 



^-amv^w . tn{t _ z) j 



solves the Riemann-Hilbert problem 12.4f> . with w replaced by W. It is the unique solution if we add 
additional requirements at the zeros of the weight: 

Y(z) = O ( | | J , as z — > afe, z e C\Ti, k = 1, . . . , m. 

Nothing hinders performing Step 1: with H defined by (|2.5|l we put T = YH, so that T becomes holomorphic 
in C \ Ti (including the infinity) and 

r+(t)=T_(t) ^ ™(t)\ teTi< 

However, in order to get rid of the oscillatory behavior of the diagonal entries of the jump matrix the lenses 
we opened in Step 2 of Sectional are no longer valid, at least because W has singularities on Ti. Since they 
are in a finite number, we can modify this step by opening lenses inside and outside Ti, but "attached" to 
the unit circle at a k s (see Figure l3~TTl . This deformation of the contour makes the definition of K in H2.6(l 
consistent (after replacing w with W), in such a way that U = TK has the jumps 



U + {t) = U^(t)Ju{t), te 7i uT 1 u 7o , 
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0_ 



fin 



9. 4 



02 



Fig. 3.1. Opening lenses. 



with 



Ju(t) 



if t G Ti, 




(3.2) 





-1/W(t) 
1 

1 

\/{t n W{t)) 

Step 3: Our next goal is to handle the jump on Ti via the global parametrix N built in Section fusing 
the Szego function. However, for W as in (|3.1|l the Szego function is, in general, no longer single-valued in 
the neighborhood of Ofc's, and a short digression is convenient in order to discuss briefly the form of this 
function and its multivaluedness. 

For the sake of brevity we define the set of singularities of the weight, A = {a\, . . . , a m }. We fix for 
what follows < S < 1 — p, such that additionally 6 < ^ min^ |a, — aj\, so that all neighborhoods Bs(a,k) 
(see definition 12.2611 s ) are disjoint. Denote also 



c fe = {z e C : \z - a k \ = 5} , k = 1, 



(3.3) 

we will use the standard 



as well as B = U^_ 1 B5(afe). Furthermore, given a subset X C C and a value a € 
notation a • X = {ax : x G X}; consistently, .4 • X = U' fe n =1 (a fc • X). 

In order to construct explicitly the Szego function for the modified weight W we introduce the generalized 
polynomial 



(3.4) 



k=l 



and select its single- valued analytic branch in C \ {yiuU^ak ■ [1, +oo)) by fixing the value of g(0). With this 
convention we can write the Szego functions for the modified weight W: 



D i (W;z) = ^-D i (w;z), D c (W;z)= Dc{w;z) 



q 2 (0)q(l/z) 



a ■ 



(3.5) 
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In particular, D[(W;z) is holomorphic in Di/ P \ (U^^afe • [1,1/ p)), D c (W;z) is holomorphic in {z 6 C : 
|«|>p}\(U£ =1 a fc -(p,l]), and 

' ' D e (W;oo) = D e (w;oo)=T>0, (3.6) 




A(W;0) A(t»;0) 

where r has been defined in l|2.7(l . Furthermore, with the orientation of the cuts toward infinity we have for 
k = 1, . . . ,m: 

[D i (W;z)] + = e - 2 ^[D i (W;z)}_, z € a fc • (1, l/p) , 
[A(^;^)] + -e- 2 ^[A(^;z)]_ , zea fe -(p,l). 

By definition H1.6(l and formulas (|3.5|) we have 

S{W;z) = D i (W;z)D e (W;z) = 

that is also analytic and single-valued in the cut annulus {p < \z\ < l/p} \ (A ■ (p, l/p)); furthermore, with 
our assumptions on S, function 

S (W- z) = { e ™' 3kS ( W > z ) ' if z e B s( a k) and arg(z) > arg(a fe ) , 
\e-™^S0^;z), if zeB«(a fc ) and arg(z) < arg(a fe ) , 

is holomorphic in Bg(ak), k = 1, . . . , m. So, we can define 

tf fe d =.S fe (W;a fe ) eT, fe = l,...,m. (3.9) 

Now let us get back to the global parametrix N(z) = N(W;z), given by formula l|2.9|l . that is well 
defined, has the same jumps on Ti as U(z), and by (|3.6|) . it exhibits the same behavior at infinity. Hence, 
U(z)N^ 1 (z) tends to / as z — * oo, and is holomorphic in C \ (7 U 7;). The jump on these curves is again 
exponentially close to identity, except in a neighborhood of the zeros of the weight W . This is a new 
feature, and we have to deal with this problem separately. 

Step 4: local analysis. Let us pick a singular point et^ G A. For the sake of brevity along this subsection 
we use the following shortcuts for the notation: a = ttfc, /3 = b = B${a,k), c = Ck (where Ck were defined 
in (|3.3|l ), and b + = {z £ b : arg(z) > arg(a)}, b~ = {z 6 & : arg(z) < arg(a)}. We also write Vlj = f fij n 6, 
where j € {i, e, 0, oo}, and analogous notation for curves: Ti = Ti n 6, etc. 

The goal is to build a matrix P such that it is holomorphic in b \ (Ti U 7i U 7 C ) , satisfies across Ti U 71 U 7 e 
the jump relation P + — P-Jjj, with Jjj given in (|3.2|l . with the same local behavior as U close to z = a, and 
matching N on c. This analysis is very technical, and we refer the reader to |llj for details, and describe 
here the main ideas in a very informal fashion. 

As a first step we reduce the problem to the one with constant jumps. Let us denote Ti = a • (1 — 5, 1) 
and r c = a ■ (1, 1 + S), oriented both from z — a to infinity (see Fig. I3.2JI . Let w 1 ^ 2 (z) and z 1 / 2 denote the 
principal holomorphic branches of these functions in b, and W 1 ^ 2 {z) = q(z)q(l/z) w 1 ^ 2 {z), with q defined in 
(13.41) . Then W 1 / 2 is holomorphic in b \ a ■ (1 — S, 1 + S), and according to l|3.7|) . 

lF V2 (z) VFl /2 (z) 
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Fig. 3.2. Local analysis in b. 



Thus, if we define 



m*) = < 



e~ ni l 3 W 1 ^{z)z n / 2 , ze (Q e u fi^) n 6~, 

e -*ip W l/2( z ) z -n/2 t zG (J} i{J o ) n b+, 

K e nif3 W^ 2 (z)z- n / 2 , z e (ftiUQo) n&", 



and set 



R(z) = P(z) A(/3; zf 3 , z e b \ (A U T c U Ti U 7i U 7c ) , 



(3.10) 



(3.11) 



we get for R the following problem: R is holomorphic in b \ (r\ U T c U Tt U 7 i U 7 C ), and satisfies the jump 
relation R+(z) — R~{z)Jr{z), with 



Jr{z) = < 



1 
-1 QJ' 
1 C,\ 

1 o\ 

e 2*ifi x J » 

o e - 7Ti/3 / 

Moreover, R has the following local behavior asz^a: 



if ze T 1; 

if z e (7i n b+) u (7 C n &-), 
if ze (7in&-)u ( 7o n6+), 
if z e Ti u r c . 



R(z) = < 



(9 
(9 



z — 


af \z 


z — 


a\P \z 


z — 


a|-" | 


z — 


ah* 3 1 



1-/3 



, if z € fio U rio< 

| , if z e ^ c u hi. 
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Consider in C \ (—00, 0) the transformation 

C = -i - \og{z/a) , (3.12) 

(we omit the explicit reference to the dependence of £ from a and n in the notation) , where we take the main 
branch of the logarithm. This is a conformal 1-1 map of b onto a neighborhood of the origin. Moreover, 
Ti is mapped onto R oriented positively, T; U T c are mapped on the imaginary axis, and we may use the 
freedom in the selection of the contours deforming them in such a way that /(71) and /(7 e ) follow the rays 
{argC = ±~ ± 7r}. After this transformation we get a Riemann-Hilbert problem on the £-plane that has 
been studied for the local analysis of the generalized Jacobi weight on the real line. We take advantage of 
the results proved therein in order to abbreviate the exposition, and refer the reader to [2H| Theorem 4.2] 
where the solution <3> ((3; £) is explicitly written in terms of the Hankel and modified Bessel functions. 

Since a left multiplication by a holomorphic function has no influence on the jumps, and taking into 
account H3.11|) . we see that matrix P can be built of the form 

P(z) = E(z)*(P;Q A(/3;z)- CT3 , (3.13) 

where E is any holomorphic function in b. An adequate selection of E is motivated by the matching 
requirement P(z)N^ 1 (z) = I + C(n _1 ) on the boundary c, and is constructed analyzing the asymptotic 
behavior of the matrix-valued function <3> (/?; £) at infinity. Let us summarize the results of this analysis: 
Proposition 3.1. Let 

where Sk(W;-) has been defined in (|3.8(l . and we take the main branch of the square root. Then matrix 
P(z) = P(a,(3;z), 

P(a, (3; z) M E(z) * (/?; C) A(# z)— , (3.15) 

with C given by (|3.12|) . satisfies: 

(i) U(z)P^ 1 (z) is holomorphic in b; 
(ii) for z € c. 



i{3 ( -T- 2 a n S(W;z) 



1 



P(z)N-\z)=I+^ 1 ^ \+0 - . (3.16) 

In particular, P(z)N^ 1 (z) — I + ©(ti -1 ) for z € c. 

Step 5: asymptotic analysis. With the notation introduced in (|3.3J) and with P(a,f3;z) defined by 
(j3~14) - (j3~TK|) let us take 

P(z) = P(a k ,f3 k ;z) for z <E B s (a k ) \ (T x U 7 U 71), k = l,...,m, 

and put 

dcf (u{z)N- 1 (z), forzeC\(BUT 1 U7 c U 7i ), 



S(z) 



U(z)P- 1 (z), for Z eB\(T 1 U 7c U 7i ). 



Matrix 5* is holomorphic in the whole plane cut along 7 U C, where 
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Fig. 3.3. Jumps of S. 



7 =(7eU 7i )\5 and C^Uf =l c k 
(see Fig. 13.31) . 5(z) — > / as z — > oo, and if we orient all c^'s clockwise, <S+(i) = S-(t)Js, with 

if z G C, 



Js(t) 



(P(z)N-\z), 

1 0' 
T 2 /(z n S(W;z)) 1 
1 -z"<S(W;z)/t 2 
1 



if z e 7c \ B, 
if z e 7i \ B. 



It is clear that the off-diagonal terms of Js on 7 ; \ B and 7c \ B decay exponentially fast. On the other 
hand, by l|3.16[) . Js(z) = I + 0(l/n) for z g C. So the conclusion is that the jump matrix J$ = I + 0(l/n) 
uniformly for z € 7 U C. Then arguments such as in |?|QJ[5] lead to the following conclusion: 
Proposition 3.2. Matrix S satisfies the following singular integral equation: 



S(z)=I + 



2ni 



(J s (t)-I)dt | 1 
t — z 2ni 



(S.(t)-I)(J s (t)-I)dt 



t- z 



where we integrate along contours 7 U C with the orientation shown in Fig. \3.3[ In particular, integrating 
counter-clockwise, 



tl ^ J<* 



t 



locally uniformly for z S C \ (7 U C). 
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Now formula ((3.16(1 and the residue theorem yield for z 6 C \ (7 U B) , 

(3 k -r- 2 a n k d k \ 
-fa J 



m 

S(z)=I+-J2 



k=l 



vr 



(3.17) 



We are ready for the asymptotic analysis of the original matrix Y (and in particular, of its entries (1, 1) and 
(2,1)). 

Unraveling our transformations we have 



Y(z) 



S(z)N(z)K- 1 (z)H- 1 (z), HzeC\B, 



S{z)P{z)K- 1 {z)H- 1 {z), ifzeB. 
We must analyze the consequences of these formulas in each domain (see Fig. 13.4(1 . 




Fig. 3.4. Domains for the asymptotic analysis. 

We will do it only for the interior domain Qq \ B, where 

w M-(-r/«V;.> D,W o 2)/T > *M-*CW- 

Hence, Y(z) = S(z)N(z), so that 

T T 

Yit(z) = ; r S\v(Z) , Y<r\ (z) = ; Soiiz) . 

V ' Di(W;z) V V ' Di(W;z) V ' 

Taking into account ((3.6(1 and ((3.17(1 . and recalling that = Y\\ we obtain 

AMQ) 1 (^Mk „ +1 L «/i ' 



«n(«) = 777777 — \ _ / . a fe + ^ 

A(W;z) 7i a fe - z 
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valid uniformly in this domains. It shows in particular that for every compact set K C Di there exists 
N = N(K) G N such that for every n > N, each has at most m — 1 zeros on K, and these zeros should 
be asymptotically close to those of the rational fraction 

/c— 1 

Evaluating 3> n (z) at z = we can obtain asymptotics for the Verblunsky coefficients a„; I leave this as an 
exercise for an interested reader. 

Remark: Orthogonal polynomials with respect to non-analytic (but smooth) and non- vanishing weights 
with finite jump discontinuities on Ti were studied in |13| . using a different (but complementary) method 
based on the d problem. There are very interesting similarities with our case: roughly speaking, both zeros 
and jump discontinuities of the weight have the same effect on the asymptotics of ^n's inside the unit disk 
Di. More precisely, according to |13| . if W > is sufficiently smooth and piecewise continuous, with jump 
discontinuities at au € Ti, k = 1, . . . , to, then with a proper normalization, & n (z) for z on compact subsets 
of Di will be asymptotically close to rational fractions of the form (|3.18|> . with the basic difference that now 
coefficients Pk stand for the magnitude of the jump of log W at a k (cf. formula (56) in ^3] and the fact that 
for the scattering function on Tj, 

S(W;e ie ) =e n ^, 

where is defined by (21) in |l."t|). These similarities have as a common ground the duality of both cases: 
for W in Ij3.1|l . the imaginary part of logVF has jumps proportional to /3fc's, while in the finite jumps 
correspond to Re (log W). 

Finally, recall that the leading coefficient n n of the orthonormal polynomial ip n is expressed in terms of 
Y by l2i(0) = — 27tk^_ 1 . This yields immediately the following asymptotic formula: 

k1_, = — I I > >T + CM — 1 | . a - x. 



This result has a consequence for the behavior of the Toeplitz determinants for W. If we define the moments 



d k = 6 z- n W(z)\dz\, 
JzeT 

then the Toeplitz determinants are 

V n (W) = det[(d^)l J=0 ] . (3.19) 
It is known (see e.g. ^0 Theorem 1.5.11]) that 

T> n (W) _ 1 



.2 



Taking into account the asymptotics of K n , l|2.8|l and l|3.6f) . we arrive at 

Theorem 3.3. Under the assumption above there exists a constant k depending on W such that 

V n {W) = x {G[2-Kw]) n n^"=i t 3 * (1 + o{\)) , n -> oo . (3.20) 

This formula is in accordance with the well known Fisher-Hartwig conjecture (see e.g. |5]), proved for this 
case (but using totally different approach and giving an expression for x) by Widom in |21|. 
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